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In this paper, the lepton flavor violating ii~V'^{V'^ = /o", 0, tj) decays are studied in the 

framework of the two Higgs doublet model(2HDM) III. We present a computation of the the 7—, Z 
penguin and box diagrams contributions, and make an analysis of their impacts. Our results show 
that, among the 7— penguins, the penguins with neutral Higgs in the loop are very larger than 
those with charged Higgs in the loop. We find that the model parameter Xt^ is tightly constrained 
at the order of 0(10""^) and the branching ratios of these decays are available at the experiment 
measure. With the high luminosity, the B factories have considerable capability to flnd these LFV 
processes. On the other hand, these processes can also provide some valuable information to future 
research and furthermore present the reliable evidence to test the 2HDM HI model. 

PACS numbers: 13.35.Dx, 12. 15. Mm, 12.60.-i 



I. INTRODUCTION 

The flavor physics is always the hot subject in particle physics. Recently, with rapid development of neutrinos 

ri 

experiment fl], the lepton flavor violation(LFV) processes of charged- lepton sector have attracted many people's 
attention. In the standard model(SM), the LFV processes are forbidden. Hence, the LFV decays are expected to 
be a powerful probe to many extensions of the SM with new LFV source and/or new particles. 

The LFV r decays have become a seeking goal in experiment. Due to the comparability of e~^e~ bb and 
e+e" T^T^ cross section (cr ~ 0.99n6) around the T(4s) energy region, large events of r leptons are available 
at BaBar and Belle {CsaBar 'i70fb^^,CBeiie 710/5^^). And now the tau pairs production has attained the reach 
of 10~^. The tau factory has performed the experimental search for the tau radiative decays and r — s- 3/ decays. 



as well as r — > decays 



2|. The current experimental upper limits of the r — > /i /O°(0, decays with 543/6 ^ 



of data at Belle laboratory are 3[: 



B{t- fi-p°) < 6.8 X 10-^ 90%CL 
B{t- fi-(j)) < 1.3 X 10"^ 90%CL 
B{t- fi-Lu) < 8.9 X 10"^ 90%CL 



(1) 



There are also lots of theoretical researches on t ^ IV'^ decays in many possible extensions of the SM. For 
example, Saha et al. have deliberated constraints on the parameters from r lp^{(f), K*'^ , K**^) decays in RPV 
SUSY model J|. Ilakovac et al. found only the ratios of e^p"(0, tt") decays reach the order of 10^^ in models 

with heavy Dirac or Majorana neutrinos 5|. The case of r ^ IP{V'^) decays in topcolor model have been considered 
by Yue Chongxing et al.[6]. Such investigations also have been presented in MSSM and minimal susysemmetry 
50(10) modelsl?!, a general unconstrained MSSM model 8| and two constrained MSSM seesaw models[9] as well. 

In our previous work 1^, we have studied the r fJ-PiP — tt", rj, t]') decays in 2HDM model III. In this model, 
there exist flavor-changing neutral currents (FCNCs) at tree level. In order to satisfy the current experiment 
constrains, the tree-level FCNCs are suppressed in low-energy experiments for the first two generation fermions. 
While processes concerning with the third generation fermions would be larger. These FCNCs with neutral Higgs 
bosons mediated may produce sizable effects to the t— /i transition. The r pP decays could yield one pseudoscalar 
meson from the vacuum state through the scalar and pseudoscalar currents. Hence, this type decay could occur at 
the tree level through the neutral Higgs bosons exchange. In this paper, we extend our discussion to the case of one 
vector meson in the hadronic final state. Different from pseudoscalar meson, the vector meson is only generated 
through vector currents and therefore receive no contributions of the neutral Higgs at tree level. So we consider 
the effects with Higgs bosons in the loop. There are the 7—, Z penguin and the box diagrams for the r — > ^"(0, uj) 
decays. For the instance of vector meson K*'^{K*^), the LFV processes could occur at loop level likewise but the 
additional loop at the hadronic vertex would generate one suppressed factor. So these two decays are not discussed 
in this paper. Our results suggest that, in the 7— penguins, the contributions of penguin with neutral Higgs bosons 
in the loop is greater than those of penguin with charged Higgs bosons in the loop. The model parameter is 
restrained at 0(10^"^) and the decay branching ratios could as large as the current upper limits of 0(10^^). For 
—t iJL^ PP processes, we will make further study in our later work. 

The paper is organized as follows: In section II, we make a brief introduction of the theoretical framework for 
the two-Higgs-doublet model HI. In section HI, we present the decay amplitudes and the numerical predictions for 
the branching ratios. Our conclusions are listed in the last section. 



II. THE TWO-HIGGS -DOUBLET MODEL III 



As the simplest extension of the SM, the Two-Higgs-Doublet Model has an additional Higgs doublet. In order 
to ensure the forbidden FCNCs at tree level, it requires either the same doublet couple to the it-type and d-type 
quarks(2HDM I) or one scalar doublet couple to the tt-type quarks and the other to d-type quarks(2HDM II). 



While in the 2HDM III 



11 



12l |. two Higgs doublets could couple to the w-type and rf-type quarks simultaneously. 
Particularly, without an ad hoc discrete symmetry exerted, this model permits flavor changing neutral currents 
occur at the tree level. 

The Yukawa Lagrangian is generally expressed as the following form: 



(2) 



where Hi{i = 1,2) are the two Higgs doublets. Qi.L is the left-handed fermion doublet, Uj^R and Dj^R are the 
right-handed singlets, respectively. These Qi.L, Uj^R and Dj^R are weak eigenstates, which can be rotated into mass 
eigenstates. While rj^'^ and S,'^'^ are the non-diagonal matrices of the Yukawa couplings. 
We can conveniently choose a suitable basis to denote Hi and H2 as: 



Hi 



1 

V2 




Ho 



1 



V2H+ 



\ 



V2 



02 + J 



(3) 



where G"'^ are the Goldstone bosons, H^ and are the physical charged-Higgs boson and CP-odd neutral Higgs 
boson, respectively. Its virtue is the first doublet Hi corresponds to the scalar doublet of the SM while the new 
Higgs fields arise from the second doublet H2. 

The CP-even neutral Higgs boson mass eigenstates H° and h° are linear combinations of 0° and 0° i'^ Eq.(l3]), 

= 05 cos a -I- 02 sin a, ft," = —05 sin a + 02 cos a, (4) 



where a is the mixing angle. 

After diagonalizing the mass matrix of the fermion fields, the Yukawa Lagrangian becomes 



i3i 



Ly = -UMuU - DMdD + -x" (UMu-f5U - DMDI5D) 

V 



^X'DV^KM WuR - MdL] U - ^x^UVcKM [MdR ~ MuL] D 



V2 



U + D 



D 



H°— \ \/2 

-^U < — M^ cos a + 
V2 [ V 

—;=U < Mu sin a - 

V2 V 



^R + ^^L 



I H°— \ \/2 

sin a > U —D < — Afn cos a 4- 

J V2 \ V 

COS a > U i=D < — Md sin a 

V2 \ V 



sin a > D 



COS a>D 



H+U 



VckmS,^ R — ^^VckmL 



D-H-D 



^^''^CKM^ ~ ^CKM^R 



U 



(5) 



where U and D now are the fermion mass eigenstates and 



where V^'j^ are the rotation matrices acting on up and down- type quarks, with left and right chiraUties respectively. 

Thus VcKM = i^L V^F ^^'^ usual Cabibbo-Kobayashi-Maskawa (CKM) matrix. In general, the matrices t)^'^ 

of Eq.® are diagonal, while the matrices are non-diagonal which could induce scalar-mediated FCNC. Seen 

from Eq.lO, the coupling of neutral Higgs bosons to the fermions could generate FCNC parts. For the arbitrariness 

of definition for ^^^'^ couplings, we can adopt the rotated couplings expressed £,^'^ in stead of hereafter. 

11 

In this work, we use the Cheng-Sher ansatz|ll| 



- 

which ensures that the FCNCs within the first two generations are naturally suppressed by small fermions masses. 
This ansatz suggests that LFV couplings involving the electron are suppressed, while LFV transitions involving 
muon and tau are much less suppressed and may lead to some loop effects which are promising to be tested by the 
future B factory experiments. In Eq.®, the parameter is complex and i,j are the generation indexes. In this 
study, we shall discuss the phcnomcnological applications of the type III 2IIDM. 

III. THE DISCUSSION FOR r ^ ^ I/" DECAYS 

As we have mentioned above, one vector meson could not be generated from the vacuum state through the scalar 
and/or pseudoscalar currents. In 2IIDM model III, the neutral Higgs bosons mediated tree and penguin diagrams 
have no contributions to /i^V^" processes. Accordingly, their decay amplitudes acquire contributions from 

the 7—, Z— penguin and box diagrams. Comparing to the r fiPdecays, in addition to neutral Higgs bosons, the 
penguin with charged Higgs bosons in the loop also contribute to these decays. We will make a detail analysis of 
their effects in the later paragraphs. The penguin diagrams at the quark level pertinent to these decays are list in 
Fig.l. 

The amplitudes could be factorized into Icptonic vertex corrections and hadronic parts described with hadronic 
matrix elements. In dealing with hadronic matrix elements, we take the generalized factorization approach and 
write the hadronic matrix elements as < ^^157^510 >= —my /ye* with the decay constant fy. The quark contents 
of p'^ meson are chosen as = '^{~'^^ + dd). For the vector (f> — lu meson system, we employ the ideal mixing 
scheme between (/)(1020) and cj(782) which is supported by existing data: (j) — —ss,uj — -^{uu + dd)[l^. Then, 




FIG. 1: The 7 and penguin dia grams for r jj. qq decay, where the neutral and charged Higgs bosons are in the loop. 



the total amplitudes could be expressed as: 



M{t- 



-^-^ ■ fl ■ [Fl^ ■L + F2^-R + F3^- jpL + FA^ ■ jpR] ■ r 



1 



1 



(^{V'l-uYu - -dYd - -sYs\0) 



Mi 



Mb 



8m^ 



F^-L + F^ ■R + Fi ■ jpL + Fl ■ 7pi?] • r 



2 



TO 



W 



F^°^ ■ L + F^" • R + i^3«°" • 7pL + F^" • jpR + Fg^"" • iapxL + Fg""" • lapxR 



^box 



pbox 



box 



ybox 



ybox 



(9) 



where M-yjMz and Mtox are the amplitudes of the 7— penguin, Z penguin and box diagrams. The relevant 
auxiliary functions are listed in Appendix. 

In our calculation, the input parameters are the Higgs masses, mixing angle a, \Xij\ and their phase an gles 6. 



Given the constraints from the current experiment permits and theoretical considerations 



fl 



15 



16 



17 



18 



18 



20|, 



we assume 



toh± = 200GeV, mno = leOGe^, to,,o = llbGeV, m^o = UOGeV, a = 7r/4, 
|A„„| = 150, |Add| = 120, |A,,| = 10, lAttl = lAtcl = |A„t| = 0.03, 

|A,,| = \Xbb\ = \Xdb\ = \Xbs\ = 100, 9 = 7r/4, (10) 



where the Higgs masses satisfy the relation 115GeV < to^o < uiao < rriffo < 200GeV 
value of Xtt ■ Xbb is approximate to three 



17 



18|. 



15 



16 



ITj l , and the absolute 



Using the above parameters, we could get the contributions of three diagrams to these decays. As we expected, 
the contributions of box diagrams are O(10~^^) order or so which are very smaller than those of 7— and Z— 
penguins. Hence, we neglect the box diagrams contributions. We have studied the relation of branching ratio and 



TABLE I: Constraints on the X-r^ from r /i p^{4>,ui) decays in the 2HDM III. 



Decay modes 


Bounds on A-r^ 


Previous Bounds 




< 1.26 X 10"^ 


~ 0(1) [21] 




< 2.45 X 10"^ - 


-O(IO) [10, 22, 2^, A^^ ~ 0(10) - O(IO^) [16] 


r )iuj 


< 1.48 X 10"^ 


A^f. ~ 0(10^) - 0(10^) [24] 



Xrfj,- The computation indicate that the variation of 0,-^, the phase angle of parameters A,-^, does almost not affect 
the values of branching ratios. So we take 0t-^ = 7r/4 as literatures do. 

The Fig. 2 gives the total penguin contributions denoted by the solid line. We denote the 7— penguin and the Z— 
penguin contributions by the dash line and the dot line, respectively. Due to the suppressed factor 0{l/ni^) from 
the Z propagator, the Z penguin contributions arc supposed to be lower than those of the 7— penguin. These decay 
amplitudes have common leptonic parts, so the differences of decay amplitudes mainly come from the hadronic 
parts. For the similar contents of p° and uj, the curves of M P*' *md — s- decays display similar trend, 

namely, their Z penguin contributions are lower one order than those of their 7— penguin. However, for ^.^ (f) 

decay, the magnitudes of Z penguin are close to the 7— penguin contributions. 

The relations of ratios versus \Xt^ji\ are also presented in Fig. 2, where the horizon lines denote the experimental 
upper limits. Evidently, one can see from Fig. 2 that these branching ratios rise with the increase of JA^^]. We have 
got the constraints on JAr^^] from the experimental data, which are list in Table. I. It is obviously that the parameter 
\\T^l\ is restrained at the order of 0(10^''). The JAr^J constraints for jj.^ p'-^uj) are little severe than that of 

T~ 4) decay. The bounds of JAr^J from different phenomenological considerations 



21 



22 



23 



2M are 



demonstrated in Tab. I, too. Comparing the values of A,-^ in Tab. I, one can see that our constraint is stringenter 
than the limits in literatures. 

Now we illustrate the contributions of 7— penguin with neutral and charged Higgs bosons in the loop. In Fig. 3, 
the solid line denotes the 7— penguin contributions, the dash line and the dot line denote the contributions of 7— 
penguin with neutral Higgs bosons in the loop and those of 7— penguin with charged Higgs bosons in the loop, 
respectively. Apparently, the contributions of 7— penguin with neutral Higgs in the loop are quite higher by nearly 
four order magnitudes than those of 7— penguin with charged Higgs in the loop. And the dot line and the solid 
line coincide with each other for three decays. As a result, the contributions of 7— penguin with neutral Higgs in 
the loop arc dominated one. 

The contributions of Z penguin with charged and neutral Higgs bosons in the loop are demonstrated in figure 
4. The solid line denotes the Z penguins contributions, the dash line and the dot line denote the contributions 




K\ K\ K\ 

(a) (b) (c) 

FIG. 2: The branching ratios versus model parameter \\th\ with Ot/j, = 7r/4, (a) for t~ — > At~p'' decay, (b) for t~ — > n~ 4) 
decay, and (c) for r~ — » decay. The solid line denotes the total contributions; the dash line and the dot line denote the 
7— and Z penguin contributions, respectively. The horizontal lines are the experimental upper limits. 




(a) (b) (c) 



FIG. 3: The branching ratios versus model parameter |At^i| with drii = 7r/4, (a) for M P*^ decay (b) for fj,~cf> 

decay, and (c) for r~ — » fj,~u! decay. The solid line denotes the 7 penguin contributions; the dash line and the dot line 
denote the contributions of 7— penguin with neutral Higgs bosons in the loop and those of 7— penguin with charged Higgs 
bosons in the loop, respectively. 

of Z penguin with charged Higgs bosons in the loop and those of Z penguin with neutral Higgs bosons in the 
loop, respectively. Unlike the case of 7 penguin, the contributions of Z penguin with neutral Higgs in the loop are 
rather smaller by nearly eight order magnitudes than those of Z penguin with charged Higgs in the loop. So the 
contributions of Z penguin with neutral Higgs in the loop are subordinate one. In a word, the 7— penguin with 
neutral Higgs in the loop plays a main role in these decays. 




0.0 2.0x10'' 4.0x10' 6.0x10'' 8.0x10' 0.0 3.0x10'' 6.0x10'' 9.0x10'' 0.0 2.0x10'' 4.0x10'' 6.0x10'' 8.0x10'' 

K\ K\ 

(a) (e) 

FIG. 4: The branching ratios versus model parameter |At^| with 9tij, = 7r/4, (a) for t~ — > n~ p" decay (b) for t~ — » fx~ 



decay, and (c) for r~ — > fj,~u! decay. The sohd hne denotes the Z penguin contributions; the dash hne and the dot Une 
denote the contributions of Z penguin with charged Higgs bosons in the loop and those of Z penguin with neutral Higgs 

bosons in the loop, respectively 

IV. CONCLUSION 



In summary, we have calculated the branching ratios of r 



M P 



decays in the model III 2HDM. 



Comparing to the t~ P decays, besides the neutral higgs bosons in the loop, an additional charged Higgs 

boson in the loop offer contributions to r~ — > ^~V^ decays. The impacts of the 7— penguin, Z— penguin and 
those of two types Higgs in loop are formulated. It is concluded that the 7— penguin with neutral Higgs bosons in 
loop arc dominated in the 7— penguin, while the Z— penguin with charged Higgs bosons in loop mainly contributes 
to the Z— penguins. Our work suggests that the parameter \\TiJt\ is constrained at the order of O(10~^). And 
in the rational parameters space, the Br{T~ — > /U~ V) can reach the experimental values. With the experiment 
luminosity increasing, these LFV decays are available to the collider's measure capability. Our study is hoped to 
supply good information for the future experiment and explore the structure of the 2HDM HI model. 

Appendix 



For simplicity, we only list the amplitude of 7— penguin. 
The amplitude of 7— penguin diagrams is 



2ml 



jl ■ [Fl^ ■L + F2^- R + F3^- jpL + FA^ ■ jpR] ■ r 



(11) 



Where the auxihary functions Fj are written as: 



F2^ 



FU = 



fc2 



.1 .l-x 

■ I dx I I 

Jo Jo 



mrXruKr ' X 



rrirX 



■ K^I-XtT ■ X 1 



Sc{x,y,mjj_,xtc) 2^ 



A;2 



1 p1-x 

dx / dy 
Jo 



Sc{x,y,mjj.,xtc) ' 2 
"^At • K^i^TT -y rrir ■ X 



Sc{x,y,m]j_,xtc) 



+ - 



™M • X%f,Xrr ■ y 



1 



Sc{x,y,m]j_,xtc) 2 



X 2/ -PW - 



X X \ -p1 



(12) 



(13) 



= 



FA^ 



rriTy/mTm^ ( 1 



fc2 



(m2 - m2 ) Jo 



fc2 



"^At • X*^Xrr{x- l)ln 



.56(:c) 



2 



dx 

^0 



+ 



(m2 - ml) Jo 

Jo Jo ^J■ ) 



Xl^Xrrix - 1) • ln5„(a;,xte) - m;in56(x) +-M, 



dyQij , 

(14) 

+N* 

(15) 



The foUowings are expressions of Ji, Ki, M^, N and Qj. 



Mho 



Mao 



- J AO — 



2iImXT 



(3;, y, TO^o , ) {x , y , rn^o , Xf^) 

Kho = {rrirX* + m^Xr^,) x J^o, Kf^o = (m^A* + m^Xr^,) x J*o, ATao = 



AO ' -^tn 

A;^(m^Ar^-m^A* ) 



[x{mlu}*Xr^ + nirmfj,ojX*i^) - UsimlXr^ + m^m^A*^)] InS'f (x,x^ ) 

-[x{m1u;*Xr^, + m^m^wA*^) - (m^w* + ml(j)Xri^ - m^m^WgA*^] InS'f °(x) 

[ dx [ 
Jo Jo 



dy{ uj Xr^[\n h ^ 



+ 



mi-m^A*^[(x + y)uj - Us] - mltoXr^ 



{xjyjTnjjojX^ ) 



{x,y,m 



HO ' tn 



H0\ 



, 1,0, 



M^o = [x{m^v*Xrix + nirm^vX* ) - Vs{m^Xr^ + m^m^A* )] InS^ {x,Xt„) 



-[x{mlv*Xru. + rrirm^vX*^) - {m^v* + mlv)Xr^, - nirmi^VsX* ] lnS'^°(x) 



Qho = [In 



Sc{x,y,m]j_,xtc) , m'^ix^ - X - 1) - mly ^ 



+ 



(x, y, m^o) arj'^) 



+ 



mrm^A;^[(x + y)v - Vs] - mlvXrn 



]w*A. 



{x , y , mf^o 1 Xf^) 

x{mlX*^Xr^, + mrm^XrrX*f,) - 2iImXrT{fnrm^X*^ - m^A^,,)] ln5^°(a;,x^°) 



Qao = [In 



[x(m2 a;^A^^ + mrm^XrTX%^) - (m^ a;^ - mlXrr)Xrix - 2i'mrm^X*i^ImXrT] lriS^°ix) 
Sf{x,y,m\o,xf°) m2(x2-x- 1) -m^y 



>5c [x, y, rrij^o, Xj„ ) 



+ 



mlXTTXrn + m^m^A;^[(x + y)Arr - 2i/mArr] 



-]a;.a. 



{x, y, fn^Ao ' "^tn ) 



oj — (Arr sin^ a + sinacosa), u = (Ai-r cos^ a — sinacosa) 
ojs — 2 sin^ aReXrr + cos 2a, Vs — 2 cos^ aReXrr — cos 2a 

The integrate function expressions are : 

771 

Sa{x,xtc) = (x - l)(xtcX - 1), Sb{x)^l-x, 5c (x, ?/, , x^c) = m^- [x + (x^ - X + y)xtc] , x^c ^ 2"" ^ 

2 

771 

'S'a(a;,a;J„) = (x - l)(xj„x - 1) £'^(2:) = ^^(a;), S'^(x, m^, a;J„) = 77i2[y + x + 2;J„x(a; - 1)], xj„ = ^ (16) 

mi 
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